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Abstract 



The possible occurrence of highly deformed configurations in the Ca 
di-nuclear system formed in the 28 Si + 12 C reaction is investigated by an- 
alyzing the spectra of emitted light charged particles. Both inclusive and 
exclusive measurements of the heavy fragments (A > 10) and their associated 
light charged particles (protons and a particles) have been made at the IReS 
Strasbourg VIVITRON Tandem facility at bombarding energies of -E; a b( 28 Si) 
= 112 MeV and 180 MeV by using the ICARE charged particle multidetec- 
tor array. The energy spectra, velocity distributions, and both in-plane and 
out-of-plane angular correlations of light charged particles are compared to 
statistical-model calculations using a consistent set of parameters with spin- 
dependent level densities. The analysis suggests the onset of large nuclear 
deformation in 40 C a at high spin. 
PACS number(s): 25.70. Gh, 25.70. Jj, 25.70.Mn, 24.60.Dr 
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I. INTRODUCTION 



The formation and binary decay of light nuclear systems in the Acn < 60 mass region 
produced by low-energy (Ej & < 7 MeV/nucleon) heavy-ion reactions has been extensively 
studied both from the experimental and the theoretical points of view []]]. In most of the 
reactions studied the binary breakup of the compound nucleus (CN) is seen as either a fusion- 
fission (FF) or a deep-inelastic (DI) orbiting || process. The large-angle orbiting yields 
are found to be particularly strong in the 28 Si + 12 C reaction M, as illustrated by Fig. 1 which 
summarizes some of the experimental results that have been collected for this system; i.e., 
orbiting cross sections [U,U and total evaporation residue (ER) cross sections HBHl^l- Since 
many of the conjectured features for orbiting yields are similar to those expected for the FF 
mechanism, it is difficult to fully discount FF as a possible explanation for the large energy- 
damped 28 Si + 12 C yields [|3|-|5|. However, FF calculations |3J significantly underpredict the 
cross sections measured in the carbon channel by almost a factor of 3, thus suggesting an 
alternative mechanism (see Fig. 1). FF, DI orbiting, and even molecular-resonance behavior 
may all be active |13[ in the large-angle yields of the 28 Si + 12 C reaction [|14],|]1|]. The back- 
angle elastic scattering of 28 Si ions from 12 C displays structured excitation functions and 
oscillatory angular distributions in agreement with the relatively weak absorption of this 
system |16|| . Moreover, the resonant gross structure jn[] is fragmented into very striking 
intermediate width resonant structure . 

Superdeformed (SD) rotational bands have been found in various mass regions (A = 
60, 80, 130, 150 and 190) and, very recently, SD bands have also been discovered in the 
N = Z nuclei 36 Ar |17|JT8|1 and 40 Ca ||19|| . These new results make the Aqn ~ 40 mass 
region of particular interest since quasimolecular resonances have also been observed in 
both the 36 Ar and 40 Ca dinuclear systems [13|]. Although there is no experimental evidence 
to link the SD bands with the higher lying rotational bands formed by known quasimolecular 
resonances, both phenomena are believed to originate from highly deformed configurations 
of these systems. Since the detection of light charged particles (LCP) is relatively simple, 
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the analysis of their spectral shapes is another good tool for exploring nuclear deformation 
and other properties of hot rotating nuclei at high angular momenta. Experimental evidence 
for angular-momentum-dependent spectral shapes has already been extensively discussed in 
the literature [p0|-p4| and, in particular, the 24 Mg + 16 reaction |[35|| , which reaches the 
40 Ca CN, has been studied in detail. Strong deformation effects have also been deduced from 
angular correlation data for the fusion reaction 28 Si( 12 C,2a) 32 S 5 . s . at Ei a b = 70 MeV |]36 



We decided to investigate the 40 Ca nucleus produced through the 28 Si + 12 C reaction at 
beam energies of Ei a b = 112 MeV and 180 MeV. As can be observed in Fig. 1, model 
calculations suggest that at the lowest incident energy of the present work the orbiting 
process is dominant for the C and O channels whereas at Ei a b = 180 MeV a large part of 
the O and N fully-damped yields may also result from a FF mechanism. In this article we 
will focus on the LCP's found in coincidence with heavy fragments. 

The present paper is organized in the following way. Sec. II describes the experimental 
procedures and the data analysis. Sec. Ill presents the inclusive and the exclusive 28 Si 
+ 12 C data (part of the experimental results presented here in detail have already been 
briefly reported elsewhere PT|-[4"T|] ) . The data are analyzed using the Hauser-Feshbach 
evaporation code C AC ARIZ O pnJ2T]J21] using a consistent set of parameters which has 



been found to successfully reproduce 24 Mg + 16 reaction results [p5|1 . The full statistical- 
model calculations, using Monte Carlo techniques to account for the experimental acceptance 
when comparing to the experimental exclusive data, are presented in Sec. IV. The strong 
cluster emission of 8 Be which is observed by the experiment is also discussed in this section. 
We end with a summary of our conclusion in Sec. V. 



II. EXPERIMENTAL PROCEDURES AND DATA ANALYSIS 

The experiments were performed at the VI VITRO N Tandem facility of the IReS Stras- 
bourg laboratory using 112 MeV and 180 MeV 28 Si beams which were incident on 12 C targets 



(160 and 180 /ig/cm 2 thick, respectively) mounted in the ICARE scattering chamber |f44] , f45 



The effective thicknesses of the 12 C targets were accurately determined using Rutherford 
back scattering (RBS) techniques with 1 H and 4 He beams provided by the Strasbourg 4 MV 
Van de Graaff accelerator p9| , ^2| , p3| . The carbon buildup corrections were found to be less 
than 2% of the total of C atoms in the targets. Both the heavy fragments (A > 10) and their 
associated LCP's (protons and a particles) were detected in coincidence using the ICARE 
charged-particle multidetector array [£|4],£|I| which consists of nearly 40 telescopes. Inclusive 
data have also been collected for heavy fragments and LCP's. 

The setup of the measurement at E; a ft( 28 Si) = 112 MeV was designed to collect only 
in-plane coincidences, whereas the setup at E/ a &( 28 Si) = 180 MeV allowed both in-plane 
and out-of-plane angular correlations to be measured (see Table I). The heavy fragments 
consisting of ER as well as quasi-elastic, deep-inelastic, and fusion-fission fragments were 
detected in 8 gas-silicon hybrid telescopes (IC), each composed of a 4.8 cm thick ionization 
chamber, with a thin Mylar entrance window, followed by a 500 /im thick Si(SB) detector. 
The IC's were located at forward angles in two distinct reaction planes (for each plane, the 
positive and negative angles in Table I are defined in a consistent manner as for the LCP 
detectors described below). The in-plane detection of coincident LCP's was done using 4 
three-element telescopes (TL3) (40 /im Si, 300 /im Si, and 2 cm CsI(Tl)) placed at forward 
angles, 16 two-element telescopes (TL2) (40 /tm Si, 2 cm CsI(Tl)) placed at forward and 
backward angles and, finally, two other IC telescopes located at the most backward angles 
(see Table I). The CsI(Tl) scintillators were coupled to photodiode readouts. The IC's were 
filled with isobutane at a pressure of 30 Torr for the backward angle telescopes and of 60 
Torr for the forward angle detectors, thus allowing for the simultaneous measurement of 
both light and heavy fragments. 

For the measurement at E; a b( 28 Si) = 180 MeV, two distinct reaction planes were defined 
as shown in Table I. One for in-plane correlations and a second one, perpendicular to the 
LCP detection plane, for out-of-plane correlation measurements. The heavy fragments were 
detected in 10 IC's located at forward angles for both the in-plane coincidences and out- 
of-plane coincidences. Both the in-plane and out-of-plane coincident LCP's were detected 
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using 3 TL3's placed at forward angles and 24 TL2's placed at forward and backward angles. 
The IC's were filled with isobutane at a pressure of 60 Torr. 

The acceptance of each telescope was defined by thick aluminum collimators. The dis- 
tances of these telescopes from the target ranged from 10.0 to 30.0 cm, and the solid angles 
varied from 1.0 msr at the most forward angles to 5.0 msr at the backward angles, according 
to the expected counting rates. 

The energy calibrations of the different telescopes of the ICARE multidetector array 
were done using radioactive 228 Th and 241 Am a-particle sources in the 5-9 MeV energy 
range, a precision pulser, and elastic scatterings of 112 MeV and 180 MeV 28 Si from 197 Au, 
28 Si, and 12 C targets in a standard manner. In addition, the 12 C( 16 0,a) 24 Mg* reaction at 
Eza& = 53 MeV [3J| was used to provide known energies of a particles feeding the 24 Mg 
excited states, thus allowing for calibration of the backward angle detectors. The proton 
calibration was achieved using scattered protons from C5H4O2 Formwar (polyvinyl formal) 
targets bombarded in reverse kinematics reactions with both 28 Si and 16 beams. On an 
event-by-event basis, corrections were applied for energy loss of heavy fragments (A > 10) 
in the targets and in the entrance window Mylar foils of the IC's and thin Al-Mylar foils 
of the Si diodes. A correction was also applied for the pulse height defect in the Si de- 
tectors. The IC energy thresholds and energy resolution for heavy fragments are better 
than 1.5 MeV/nucleon and 0.7%, respectively. The total energy resolution of 8.78 MeV a 
particles from thorium sources has been found to be better than 2.2% for both the three- 
element and two-element light-ion CsI(Tl) telescopes. Absolute cross sections of inclusive 
measurements could be obtained within 10-12% uncertainties, with 3-5% uncertainties in the 
target thickness and to 8-10% uncertainties in the electronic deadtime corrections. More 
details on the experimental setup of ICARE and on the analysis procedures can be found 



in Refs. P2|,^3|,|38|-4l0|1 and references therein. 



III. EXPERIMENTAL RESULTS 
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A. Inclusive data 



The fragments with Z = 15-17 have inclusive energy spectra that are typical of ER's. 
The inclusive data of the binary fragments with Z = 5-14 were obtained by integrating their 
l/sin# c . m . angular distributions. In particular the results of the cross sections with Z = 
6-8 are compared with the previously measured excitation functions in Fig. 1. The 
agreement between both sets of data is satisfactory within the error bars. 

The LCP inclusive data can shed some light on the reaction mechanism. Typical inclusive 
energy spectra of a particles are shown for the 28 Si + 12 C reaction at Ei a b = 112 MeV 
in Figs. 2-(a) and 2-(b) and at E/ a {, = 180 MeV in Figs. 2-(c) and 2-(d) at the indicated 
angles. The solid points (with error bars visible when greater than the size of the points) are 
the experimental data whereas the solid and dashed lines are statistical model calculations 
discussed in Sec. IV.B. The a-particle energy spectra have a Maxwellian shape with an 
exponential fall-off at high energy which reflects a relatively high effective temperature 
(T ' s i ope « [8E* CN /A CN } 1/2 = 3.67 MeV for E lab = 180 MeV) of the decaying nucleus. The 
spectral shape and high-energy slopes are also found to be essentially independent of angle 
in the cm. system. These observations suggest a statistical de-excitation process arising 
from a thermalized source such as the 40 Ca CN and are consistent with a previous study at 
E lab = 150 MeV g§. 

The velocity contour maps of the LCP Galilean-invariant differential cross sections 
(d 2 (r/df2dE)p _1 c _1 as a function of the LCP velocity are known to provide an overall pic- 
ture of the reaction pattern. From this pattern the velocity of the emission source can be 
determined in order to better characterize the nature of the reaction mechanism. Such typ- 
ical plots [^,^,0,0 (not shown here) can again be understood by assuming a sequential 
evaporative process and successive emission sources starting with the thermally equilibrated 
40 Ca* CN and ending with the final source characterized by a complete freeze-out of the 
residual nucleus. The Galilean-invariant cross section contours form arcs that are centered 



at V c m . [HH], as expected for a complete fusion-evaporation (CF) mechanism followed by 
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isotropic evaporation. 



B. Exclusive data 

1. LCP energy spectra 

The exclusive LCP events are also largely consistent with a CF mechanism with the 
light particles being in coincidence with ER's with Z = 15-17. The data taken with the 
IC's located at more backward angles (larger than -15°) are not considered in the following 
analysis since the statistics for fusion-evaporation events are too low. The experimental 
data of Figs. 3 and 4 are given by the solid points, with error bars visible when greater 
than the size of the points. At E^ = 180 MeV, the spectral shapes of a particles in 
coincidence with P (Z=15) ER's, shown in Fig. 4, are very similar to the inclusive energy 
spectra of Fig. 2. On the other hand the energy spectra of a particles in coincidence with 
S (Z=16) ER's are more complicated as they show other non-evaporative sub-structures 
(their angular dependence are indicative of a binary decay origin) which are superimposed 
on the "statistical" Maxwellian shape. These additional non-statistical components will be 
discussed as arising from a 8 Be cluster emission in Sec. IV. C. At E iab = 112 MeV, shown 
in Fig. 3, for Z = 15 the high-energy components showing up at the large angles may arise 
from a, 3p evaporation cascades. 

Fig. 5 presents the corresponding 28 Si + 12 C exclusive energy spectra of protons emitted 
in coincidence with individual ER's (Z = 15 and 16) at E iab = 180 MeV. Their spectral shapes 
are also Maxwellian with the typical exponential fall-off at high energy, characteristic of a 
statistical CN decay process. 

2. In-plane angular correlations 

The in-plane angular correlations of a particles and protons in coincidence with ER's 
produced in the 28 Si + 12 C reaction, are shown in Figs. 6 and 7 at E iab =112 MeV and 180 
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MeV, respectively. The angular correlations are peaked strongly on the opposite side of the 
beam direction from the ER detectors which were located at QfJ^ = -15° or QfJ^ = -10° for 
the two energies, respectively. The observed peaking of the LCP yields on the opposite side 
of the beam from the IC is the result of momentum conservation. The angular correlations 
of both the protons and a particles show the same behavior for all ER's. 



3. Out-of-plane angular correlations 

Fig. 8 displays the out-of-plane angular correlations of a particles (circles) and protons 
(triangles) in coincidence with individual ER's detected at Qi a b = 10°, produced in the 
28 Si(180 MeV) + 12 C reaction. The angular distributions have a behavior following a typical 
exp(-asm 2 (0[ ab )) shape |39| , [47| - ^9|1 , with possibly two components visible for a particles in 



coincidence with Z = 14, 15 and 16, plus a broadening of the angular correlations at backward 
angles. This broadening effect may result from a particles being able to be emitted at the 
beginning or at the end of the decay chain, where the angular momentum becomes smaller 
towards the end of the chain. As the protons cannot remove as much angular momentum as 
do the a particles the broadening effect is less significant in the proton angular correlation. 
The solid lines shown in the figure are the results of statistical-model predictions for CF 
and equilibrium decay using the Monte Carlo evaporation code CACARIZO [p0|,^li p4|, as 
discussed in the next Section. 



IV. DISCUSSION 



A. Statistical-model calculations 



The evaporation of light particles from a highly excited CN is a well known decay process 
up to very high excitation energies and spins [p0|-|5^|. The interpretation of LCP data 
requires a careful treatment of the light particle emission properties in the statistical-model 
description. Most of the available statistical-model computer codes |T|J^,^0H52| are based 
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on the Hauser-Feshbach formalism and are able to follow the CN decay by a cascade of 
evaporated LCP's and neutrons. In particular, a detailed analysis of the exclusive data 
can be undertaken by the use of Monte Carlo versions |HJ of some of these statistical- 
model codes in which the filtering of the events can reproduce the experimental conditions. 
The statistical-model analysis of the present data has been performed using the Hauser- 
Feshbach evaporation code C AC ARIZ O pOf. C AC ARIZ O is a Monte Carlo version of 



the statistical-model code CASCADE [p^l , which has evolved with many modifications 



and extensions p0 j|21| , P4| from the original code. In this program the effective experimental 
geometry of the ICARE detectors is properly taken into account. It is assumed that a single 
CN is created with a well defined excitation energy and angular momentum distribution, and 
the de-excitation chain is followed step by step and recorded as an event file. The generated 
events are then analyzed using a subsequent filtering code ANALYSIS [[21]] in which the 
locations and the solid angles of all the ICARE telescopes are explicitly specified. This 
program allows the determination of the different types of events of interest. Such events 
can be sorted (singles events, coincidence events, etc.) and the corresponding particle spectra 
and angular distributions can be created. 

The CN angular momentum distributions needed as the primary input for the calcula- 
tions were specified using the critical angular momentum for fusion \j cr it and the diffuseness 
parameter AL. They were taken from ER cross section data compiled for the 28 Si + 12 C 
fusion process by Vineyard et al. [11], without including fission competition. A fixed value 



of AL = lh (optimized at low energy by a previous statistical-model analysis of this reac- 



tion ||36|| ) was assumed for the calculations. It has been checked that the calculated spectra 
are not sensitive to slight changes in the critical angular momentum or to explicit inclusion 
of the fission competition. The parameter sets used for the calculations are summarized in 
Table II. 

The other standard ingredients for statistical-model calculations are the formulations of 
the nuclear level densities and of the barrier transmission probabilities. The transmission 
coefficients were derived from Optical Model (OM) calculations using potential parameters of 
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light particle induced reactions deduced by Wilmore and Hodgson fl54f , Perey and Perey |55 



and Huizenga [|56] for the neutrons, protons and a particles, respectively. For spin regions 
where the standard rotating liquid drop model (RLDM) JST] as well as the finite-range liquid 
drop model (FRLDM) |58| still predict essentially spherical shapes, these sets of transmission 
coefficients have been found adequate in the considered mass region. However, in recent 
years it has been observed that when the angular momentum is increased to values for 
which FRLDM predicts significant deformations, statistical-model calculations using such 
standard parameters cannot always predict satisfactorily the shape of the evaporated a- 
particle energy spectra p]| -|34||. The calculated average energies of the a particles are found 
to be much higher than the corresponding experimental results. Several attempts have been 
made to explain this anomaly either by changing the emission barrier or by using spin- 
dependent level densities. Adjusting the emission barriers and corresponding transmission 
probabilities affects the lower-energy part of the calculated evaporation spectra. On the 
other hand the high-energy part of the spectra depends crucially on the available phase 
space obtained from the level densities at high spin. In hot rotating nuclei formed in heavy- 
ion reactions, the energy level density at higher angular momentum is spin dependent. 
The level density, p(E, J), for a given angular momentum J and energy E is given by the 
well known Fermi gas expression with equidistant single-particle levels and a constant level 
density parameter a: 

P(E , J) = (^ V^ ) 3/ 2 _J__ exp(2[0(£ _ A _ T _ Ej)] ^ (1) 



where T is the "nuclear" temperature and A is the pairing correction []59 |. The quantity 
Ej = 2 j 2 J(J+1) is the rotational energy, with J e jf = J x (1 + 5i J 2 + 8 2 J 4 ) being the 
effective moment of inertia, where j7o at high excitation energy and high angular momentum 
is considered to be the rigid body moment of inertia and 8\ and 82 are the "deformability 



parameters" [Eq-E2L23 



The level density parameter is constant and is set equal to a = A/8 MeV 1 , a value which 



is in agreement with previous works |35|j60|j6T[| . In principle, the value of a may be affected by 
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dynamical deformation: rotation induces rearrangement of the single-particle level scheme 



and the altered nuclear surface area |27| affects the macroscopic energy of the system. The a 



parameter becomes more important when the nuclear deformation increases ||61|| . However, 
in the present work we assume a constant value and rather introduce deformation effects 
through the deformability parameters. A constant value of a = A/8 is in agreement with 
various authors |2^j3^|, as well as with theoretical studies by Shlomo and Natowitz []60| |, by 



Toke and Swiatecki ||61|| , and with experimental results obtained very recently in the Acn 



60 mass region |j2 



No attempt was made to modify the transmission coefficients since it has been shown 



that the effective barrier heights are fairly insensitive to the nuclear deformation |[27|| . On 
the other hand, by changing the deformability parameters 5± and 82 one can simulate the 
spin-dependent level density [^2| , p4| , p7| , p9| associated with a larger nuclear deformation, and 
thus better reproduce the experimental data. 

B. Deformation effects in 40 Ca 
In the present analysis, following the procedure proposed by Huizenga et al. |27|], we 



empirically modify the phase space open to statistical decay by lowering the Yrast line 
with adjustment of the deformability parameters so as to fit the available experimental 
data p2| , p4| . We may also take into account the fact that the deformation should be atten- 
uated during the subsequent emission processes: i.e., there is a readjustment of shape of the 
nascent final nucleus and a change of collective to intrinsic excitation during the particle- 
evaporation process. Such an analysis was suggested earlier by Blann and Komoto |M| , but 



with the assumption that the deformation is a frozen degree of freedom through the decay 
chain. Dynamical effects related to the shape relaxation during the de-excitation process 



have been incorporated into statistical-model codes p5|j28|1 . For the CACARIZO calculations 
done here, it is assumed that memory of formation details are lost after each step, with only 
the conserved quantities such as total energy and spin preserved during the decay sequence. 
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The CACARIZO calculations have been performed using two sets of input parameters: 
the first one with standard liquid drop parameters (parameter set A), consistent with the 
deformation of RLDM |57J and of FRLDM with finite-range corrections of Sierk pSfl , and 



the second one with larger values for the deformability parameters (param- 



eter set B) which are listed in Table |m[ In the FRLDM [58| the CN can be considered 
as spherical or slightly deformed at low bombarding energy, becoming strongly deformed as 
the spin increases. 

The dashed lines in Fig. 2 show the predictions of CACARIZO for 28 Si + 12 C at both 



bombarding energies using the parameter set A consistent with FRLDM deformation |58 
It can be observed that the average energies of the measured a-particle inclusive spectra 
are lower than those predicted by these statistical-model calculations. The solid lines of 
Fig. 2 show the predictions of CACARIZO using the same increased values of the spin 
deformation parameters (see parameter set B given in Table II) for both energies, and the 
agreement is considerably improved. 

The exclusive energy spectra of the a particles in coincidence with individual ER's (Z 
= 15 and Z = 16) are displayed in Figs. 3 and 4 for the two bombarding energies E; a fe = 
112 MeV and 180 MeV, respectively. It can be observed that the spectra in coincidence 
with the P residues are well reproduced by using the deformation effects [ 42jp7 ,^8 , 4l| . The 



solid lines in Figs. 3 and 4 show the predictions of CACARIZO using the parameter set 
B with S\ = 2.5 x 1CT 4 and 82 = 5.0 x 10~ 7 chosen to reproduce the data consistently 
at the two bombarding energies. On the other hand, by using the standard liquid drop 
deformability parameter set A with no extra deformation (i.e. with small values of 5i and 
62), the observed average energies from the exclusive a-particles spectra £1X6 , cLS found for 
the inclusive data, lower than those predicted |39[ by the statistical model. In this case the 



CACARIZO parameters are similar to the standard parameters used in a previous study 
of the 130 MeV 16 + 24 Mg reaction [|35|, with the use of the angular momentum dependent 
level densities. 

On inspection of Figs. 3 and 4 a large difference can be noticed when comparing the 
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calculated and experimental energy spectra associated with S residues and those associated 
with P ER's ||38|| . The latter are reasonably well reproduced by the C AC ARIZ O calcula- 
tions, whereas the model does not predict the shape of the spectra obtained in coincidence 
with S residues at backward angles (9 a > 40° at E iab = 112 MeV and 6 a > 70° at E iab = 180 
MeV). An additional, non-statistical, a-particle component is suggested in Section IV. C to 
arise from a 8 Be emission process. This is consistent with the discrepancies also observed 
at backward angles in the in-plane angular correlations of Fig. 6. For both a particles and 
protons the calculations significantly underpredict the yields at the negative angles at Ei ab 
= 112 MeV. However the agreement is more satisfactory for protons at Ei ab = 180 MeV as 
shown in Fig. 7. 

As shown in Fig. 5 CACARIZO calculations are also able to reproduce the shape of 
exclusive proton spectra for the 28 Si(180 MeV) + 12 C reaction. Compared to the a particles, 
it may be mentioned that the energy spectra of the protons do not shift as significantly as 
the spin-dependent parametrization of the moment of inertia is introduced. The statistical- 
model results using the two different parameter sets reproduce equally well the experimental 
velocity spectra and angular correlations. The statistical-model calculations displayed for 
protons in Fig. 5 have been performed with parameter set B (solid lines) including the 
deformation effects (calculated curves with parameter set A are indistinguishable from the 
solid lines). 

In order to better determine the magnitude of the influence of deformation effects in the 
CN and the residual nuclei which are suggested by our choice of statistical-model approach, 
we have proposed a very simple procedure |43| , |47| -^9[] . The effective moment of inertia is 
expressed as J e ff = |MR 2 = |M(b 2 +a 2 ) with the volume conservation condition: V = 
|-7rabc, where b and a are the major and minor axis, and c is the rotational axis of the CN. 
In the case of an oblate shape a = b and J e ff = |Ma 2 and V = |7ra 2 c. The axis ratio 
is equal to 5 = a/c = (l+c>iJ 2 -|-c>2J 4 ) 3/ ' 2 . In the case of a prolate shape a = c and Jo = 
|M(b 2 +a 2 ) and V = |yra 2 b. We obtain the equation: l + (3-7)x+3x 2 +x 3 = with x = (|) 2 
= 5 2 and 7 = 8(l+<5iJ 2 +<5 2 J 4 ) 3 . The quadrupole deformation parameter j3 is equal to j3 = 
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The effects of the Yrast line lowering (increase of the level density) due to the nuclear 
deformation and the variation of the deformation parameter (3 can be quantitatively dis- 
cussed using the values in Table [TIT] for several reactions. The values of the minor to major 
axis ratio b/a and of the deformation parameter (3 have been extracted (with 10% error 
bars) from the fitted deformability parameters by assuming a symmetric prolate shape with 



sharp surfaces [27|. The assumption of oblate shapes yields similar results within 5%. It 
is interesting to note that the deformation found necessary to reproduce the 28 Si + 12 C 
reaction results is smaller than the deformation introduced by the deformability parameter 



used by Kildir et al. ^6|, who also change the transmission coefficients. 

The solid lines shown in Figures 6 and 7 are the results of CACARIZO calculations 
for the exclusive measurements. It can be observed in Fig. 6 that for = 112 MeV the 
experimental angular correlations are well reproduced by the evaporation calculations for the 
data at the opposite side from the ER detector, and this is true for correlations with both S 
and P ER's. However the calculations fail to predict the experimental data at the same side 
as the ER detector. The question of these large yields measured at negative angles remains 
open f43j. Similarly the CACARIZO calculations reproduce in Fig. 7 the in-plane angular 
correlations of a particles (circles) and protons in coincidence with all ER's, at E/ a b = 180 
MeV, for the data on the opposite side from the ER detector. They are also able to describe 
the in-plane angular correlations of protons in coincidence with individual Z = 14 and Z = 
15 on both sides of the beam. However, the excess of yields observed at backward angles 
(0/a6 = +50° to +90°) for a particles in coincidence with S may indicate the occurrence of 
a non-evaporative process, possibly of a binary nature. The excess of yields is even stronger 
for a particles in coincidence with Si. Here the hypothesis of a non-evaporative component 
arising from the 12 C breakup may be advanced. 

The solid lines shown in Fig. 8 for the out-of-plane angular correlations are the re- 
sults of CACARIZO statistical-model predictions. Once again it can be observed that the 
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statistical-model calculations are able to reproduce the proton coincidences well, but they 
fail to describe the a-Cl coincidences and the large yields found in coincidence with P and 
Si ER's in the most forward direction. The poor reproduction of the a-particle experimen- 
tal anisotropy factor is not well understood. The angular momentum dependence has been 
tested by performing calculations with two different angular-momentum windows : 10-20/i 
and 20-30/i. Whereas for protons the anisotropy is almost constant with the L-window, for 
the a particles the anisotropy is strongly depending of the chosen L-window. Nevertheless 
the flat behavior shown around 0° is present for the two particle species. The observed 
discrepancy suggests that the assumed angular distribution of the LP, which is handled 
semi-classically , may not be adequate to describe the out-of-plane data. 

Overall, we conclude that the evaporated a particles from 40 Ca CN emission reflect signif- 
icant deformation effects. The deduced deformation is comparable to that found previously 
in the analysis of 28 Si( 12 C,2a) 32 S 9 . s . angular correlation data [|36|. The extent to which these 
effects can be reasonably well quantified is dependent on the experimental coverage and, in 
particular, on the power of the coincidence trigger. It is of particular interest to note that 
the value of (5 ~ 0.5 found for the quadrupole deformation parameter of 28 Si + 12 C (see Table 
III) might be connected with the recent observation of SD bands in the doubly-magic 40 Ca 



nucleus by standard 7-ray spectroscopy methods fll9| . Correlating large prolate deforma- 
tions in the hot CN with the presence of SD bands in 40 Ca is obviously not straightforward, 
since the deformation deduced from the LCP data averages over CN configurations, while 
the SD bands are based on one of these configurations. We made the same discussion with 
the possible comparison between LCP results for the 28 Si + 28 Si reaction [|3| and 7-ray data 
displaying very deformed bands in the doubly magic 56 Ni nucleus fl6"5 |. It is interesting to 
add that the macroscopic deformation energy of 40 Ca recently calculated within a Gener- 
alized Liquid Drop Model |Mj with shell effects (using the Strutinsky method) generates a 
second highly deformed minimum where SD and highly deformed states may survive. 
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C. Non-statistical 8 Be cluster emission 

It has been shown from the analysis of Figs. 3 and 4 with CACARIZO that additional 
non- statistical components appear to be significant at both = 112 MeV and 180 MeV 
bombarding energies. However no evidence was found for additional processes at the lower 
bombarding energies of Ei ab = 70 MeV |Q and 87 MeV |)7|. To better understand the 
origin of these components, a particle energies are plotted in Fig. 9 against the energies of 
the S residues detected at Os = -10° for the 28 Si(180 MeV) + 12 C reaction for a number 
of Qf-particle emission angles. With increasing a-particle angles an increase of the energy of 
S residues and a decrease of the a energy is observed which is consistent with kinematics. 
At a = +40°, +45°, and +50° the bulk of events in Fig. 9 are of a statistical origin, and 
consistent with CACARIZO calculations, as demonstrated in Fig. 10-(b) (for Q a = +40°). 
Another statistical-model code PACE 2 |)8[] gives similar predictions. The calculations sug- 



gest that these a particles result from a cascade of a single a, two protons, and x neutrons 
rather than a 2-a,xn evaporation process. For larger angles, the two branches, corresponding 
to the contours labeled 1 and 2, although lying outside the "statistical evaporation region", 
still correspond to an evaporation process as shown by the CACARIZO calculations dis- 
played in Fig. 10-(b) for Q a = +40° and in Fig. 10-(d) for Q a = +70°. These two branches 
1 and 2 correspond to a 2-a fusion-evaporation channel with both the a particles emitted 
respectively at backward and forward angles in the center of mass. However, at more back- 
ward angles other additional contributions, corresponding to the strong peak in the contour 
labeled 2 and in the contours labeled 3 and 4, appear with increasing significance, as shown, 
for instance in Fig. 10-(c) for Q = +70°. The corresponding "folding angles" are compatible 
with the two-body kinematics required for the 32 S + 8 Be binary exit-channel. In contrast, 
the energy correlations for the a particles in coincidence with CI and P residues (not shown) 
do not exhibit similar two-body branches, the "statistical evaporation region" is consistent 
with the CACARIZO predictions, for all the measured angles. 

Although in principle the identification of the 8 Be cluster requires the coincident detec- 
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tion and mass identification of both decaying a particles fl69| , a kinematic reconstruction 
assuming a two-body 32 S + 8 Be* process is instructive. Assuming the three-body kinemat- 
ics of a a+a+ 32 S final state, it is possible to reconstruct the momentum of the "missing" 
a-particle and, hence, to deduce the excitation energy of the intermediate 8 Be fragment. In 
Figs. 11 the deduced excitation energy spectra in this channel are presented for the con- 
tributions labeled 2, 3, and 4 in Fig. 10-(c) at the indicated Q angles. From Q a = 70° 
to G a = 85° the strongest peak appears with a very narrow width. This large component, 
which corresponds to the contribution of the contour 4 visible in Fig. 10-(c), is centered at 
the energy of the ground state of 8 Be (the relative energy of the two a particles of the 8 Be 
breaking up in flight is 92 keV) and displayed as the squared part of Fig. 11- (a). From B Q = 
55° to © a = 95 ° the main bulk of the yields from contours 2 and 3 is centered at around E* 
= 3.1 MeV with an experimental width of approximately 1.5 MeV, which values correspond 
well to the known energy (E* = 3.04 MeV) and width (r = 1500 keV) of the first 2 + excited 
level of 8 Be [TOf. The short-lived 8 Be 4+ excited level at E* = 11.4 MeV 0] is not clearly 



observed due to its very broad width (r = 3.7 MeV) and the significant a-statistical back- 
ground arising from the contribution of the contour 2. For the same reasons it is hazardous 
to assign the bumps around 15 MeV to the known 2 + doublet J7^] at E* = 16.6 and 16.9 
MeV. 

Figs. 11- (b) and 11- (c) display the reconstructed excitation energy spectra of the S binary 
fragments measured at 0^ = -10° in coincidence with a particles detected at the indicated 
Q a angles by gating either on the ground state (g.s.) contour 4 (upper panel) or the 2 + state 
contours 2 and 3 (lower panel). We have performed fusion-fission calculations (not shown), 
using the Extended Hauser-Feshbach Method 0. They fail to reproduce both the excitation 
energies of the 32 S fragments, and the yields from the contributions 2, 3 and 4 j37|. These 
contributions might result from a faster binary process governed by the a-transfer reaction 
mechanism 28 Si + 12 C — > 32 S* + 8 Be, as proposed by Morgenstern et al. W%. This conclusion 



is in agreement with previous inclusive results published in Ref. |12[|. In the cluster-transfer 



picture JF^] the reaction is characterized by a "Q-value" window centered at the so-called 



"Q-optimum" , which value can be estimated semi-classically by Q op4 = (Z 3 Z4/ZiZ2-l)Ef m -, 
where the indices 1,2 and 3,4 indicate the entrance (i) and exit channel, respectively. The 
corresponding excitation energy E* = Q 99 - Q op t, where Q gg is the ground-state Q-value 
of the reaction. In this case the expected excitation energy in the 32 S nuclei is equal to 
12.9 MeV. Fig. 11 represents the calculated excitation energy of 32 S in coincidence with 
the g.s. (b) for individual angles and by adding individuals angles (lower spectrum labelled 
Total), and with the first 2 + (E x = 3.04 MeV) excited state (c) of 8 Be, respectively. The 
strong shifts of the energy distributions can be explained by the bias effects induced by the 
kinematic coincidence acceptances. The dashed lines correspond to E* = 12.9 MeV, the 
energy expected for a-transfer reaction mechanisms. In both cases the excitation energies 
(total spectra for the coincidence with the g.s.) of 32 S are consistent with these values. In 
the same way we can also have a 8 Be-transfer reaction mechanism [0 28 Si + 12 C — > 36 Ar* + 
a. In this case the 36 Ar* ejectile has sufficient excitation energy to emit either one proton or 
one a particle. This type of "transferlike" reaction can explain the disagreement observed 
in Fig. 6 between data and CACARIZO calculation for the in-plane angular correlation 
between a particles and CI residues. 



V. CONCLUSIONS 

The possible occurrence of highly deformed configurations in the 40 Ca dinuclear system 
has been investigated by using the ICARE charged-particle multidetector array at the 
VIVITRON Tandem facility of the IReS Strasbourg. The properties of the emitted LCP's in 
the 28 Si + 12 C reaction have been analyzed at two bombarding energies E; a t = 112 MeV and 
180 MeV, and compared with a statistical model that was adopted to calculate evaporation 
spectra and angular distributions for deformed nuclei. A Monte Carlo technique has been 
employed in the framework of the well documented Hauser-Feshbach code CACARIZO. 
The measured observables such as velocity distributions, energy spectra, in-plane and out-of- 
plane angular correlations are all reasonably well described by the Monte Carlo calculations 
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which include spin-dependent level densities. The magnitude of the adjustments in the 
Yrast line suggests deformations at high spins in the 40 Ca dinuclear system that are far 
in excess of those predicted by the FRLDM. The deduced deformations are comparable to 
recent 7-ray spectroscopy data for the 40 Ca nucleus at much lower spins [|19|]. The suggested 
40 Ca shapes are consistent with predictions of the Generalized Liquid Drop Model pH 
which predicts a second highly deformed minimum in this system resulting from shell effects 
where SD and highly deformed states may survive. A component is found in the a-particle 
energy spectra measured in coincidence with S residues that is attributed to the decay of 
unbound 8 Be nuclei, although this study does not clearly establish the mechanism resulting 
in these yields. In general, to fully explore the influence of nuclear deformation on the 
reaction mechanisms and underlying nuclear structure in the mass A^40 region will require 



sophisticated particle-7 experiments (see Refs. [|73| -[76[| for instance) using EUROBALL IV 
and/or GAMMASPHERE. These studies are necessary to better understand how the large 
nuclear deformations that are apparent in the fusion studies are related to the superdeformed 
bands discovered and/or predicted in this mass region |T7|-|r9||. 
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TABLES 



ICARE setup 
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Si(112 MeV) + 12 C 



28 



Si(180 MeV) + 12 C 



Type of detector 




<K°) 


Type of detector 


en 


<f>(°) 


IC 


±15 





IC 


±10 







-20 


±20 




±10 


90 




±25 







±15 







-30 


±20 




±20 







-35 







±20 


90 




-40 


±20 




±25 







±130 













-150 











TL3 


±15 


±20 


TL3 


±30 







±25 


±20 




±35 







±35 


±20 










±45 


±20 








TL2 


16 telescopes each 5° from 


TL2 


24 telescopes each 5° from 




0=40° to 115° and 0=0° 




9= 


40° to 95° and 0=0° 



TABLE I. Experimental setup of ICARE chosen for the 28 Si ± 12 C reaction at £Jj j, 
MeV and 180 MeV. 
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Angular-momentum distribution in CN 

Critical angular momenta L cr = 21 (Ei ab = 112 MeV) and 27h (Ei ab = 180 MeV). 
Diffuseness parameter AL = l.Oh. 



OM potentials of the emitted LCP and neutrons 



(1) Transmission coefficients as denned in the text 

(2) Multiply factor of the OM radius: RFACT = 1 



Level-density parameters at low excitation: (E* < 10 MeV) 



(1) Fermi-gas level-density formula with empirical level-density parameters from Dilg et al. [59 

(2) Effective moment of inertia 9 = IFACT frigid 

with IFACT = 1. 



Level-density parameters at high excitation: (E* > 20 MeV) 

'1) Fermi-gas level-density formula with parameters from RLDM (Myers and Swiatecki |6 
'2) Level-density parameter: a = A/8 MeV -1 



Yrast line 

Parameter set A: FRLDM (Sierk (5f 



Parameter set B: 9 = 9 sp/iei . e (l + h J 2 + S 2 J 4 ) with 6i = 2.5 1CT 4 et 5 2 = 5.0 10 -7 



7-ray width (in Weisskopf units) 



(1) El: B(E1) = 0.001 

(2) Ml: B(M1) = 0.01 

(3) E2: B(E2) = 5.0 



TABLE II. Parameter sets used in the CACARIZO calculations for the 28 Si + 12 C reaction at 
Ei ab = H2 and 180 MeV. 
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Reaction 


C.N. 






5 2 


b/a 





Reference 


28 Si + 12 C 


40 Ca 


21 


2.5-10" 4 


5.0-KT 7 


1.3 


0.47 


This work 


28 Si + 12 C 


40 Ca 


26 


6.5-10" 4 


3.3-10" 7 


2.0 


0.53 


1 


16 




28 Si + 12 C 


40 Ca 


27 


2.5-10" 4 


5.0-10- 7 


1.8 


0.51 


This work 


28 Si + 27 Al 


55 Co 


42 


1.8-10" 4 


1.8-10" 7 


1.3 


0.46 


1 


29 




28 Si + 28 Si 


56 Ni 


34 


1.2-10" 4 


1.1-10" 7 


1.6 


0.49 


1 


13 




28c: , 28c; 


56 IVi 
1\ 1 


37 

O 1 


i 9.1 n — 4 

1 ,L' 1U 


1 . 1 • 1U 


1 7 


u.ou 


1 


39 




30 Si + 30 Si 




34 


1.2-10" 4 


1.1-KT 7 


1.7 


0.50 


1 


13 


1 


35 C1 + 24 Mg 


59 Cu 


37 


1.1-HT 4 


1.3-10~ 7 


1.7 


0.51 


1 


19 




32 S + 27 Al 


59 Cu 


42 


1.3-10" 4 


1.2-10" 7 


2.0 


0.53 


I 


27 




16 o + 54 Fe 


70 Se 


34 


2.5-10" 5 


3.0-10" 8 


1.3 


0.46 


1 


32 





TABLE III. Typical quantities of the evaporation calculations performed using the statisti- 
cal-model code CACARIZO as discussed in the text. The deformahility parameters are taken 
either from the parameter set B (see Table II) for 28 Si + 12 C or from similar fitting procedures 
for the other systems studied in the literature. 
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FIGURES 



Figure 1: Experimental C (solid squares), N (solid triangles), and (solid circles) cross 
sections measured in the 28 Si + 12 C reaction |||| as compared to the calculations (dotted 
curves) performed with the equilibrium model of orbiting ||. The solid curves are the 
predictions of the transition-state model |IJ . The open squares, triangles and circles are the 
present data of the C, N, and O fully-damped yields with error bars smaller than the size 
of the symbols. The full diamonds correspond to ER cross sections quoted in Refs. [B-[l2| . 



Figure 2: Inclusive energy spectra of a particles measured in the Si + C reaction 
at E lab = 112 MeV (a) and (b), and 180 MeV (c) and (d) at 6f^ p = 40° and 45°. The 
experimental data are shown by the solid points with error bars visible when greater than 
the size of the points. The solid and dashed lines are statistical-model calculations discussed 
in the text. 



Figure 3: Exclusive energy spectra of a particles emitted at the angles +40° < 
@£? P < +65°, in coincidence with individual P and S ER's detected at -15° in the 28 Si + 12 C 
reaction at E; a t =112 MeV. The experimental data are given by the solid points with error 
bars visible when greater than the size of the points. The solid lines are statistical-model 
calculations discussed in the text. 



Figure 4 Exclusive energy spectra of a particles emitted at the angles +40° < 
@fab P < + 95 °> in coincidence with individual P and S ER's detected at -10° in the 28 Si + 12 C 
reaction at E; a t = 180 MeV. The experimental data are given by the solid points with error 
bars visible when greater than the size of the points. The solid lines are statistical-model 
calculations discussed in the text. 
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Figure 5: Exclusive energy spectra of protons emitted at the angles +40° < 6^ p <+70°, 
in coincidence with individual P and S ER's detected at -10°, at the indicated laboratory 
angles, in the 28 Si (180 MeV) + 12 C reaction. The solid lines are statistical-model calculations 
discussed in the text. 

Figure 6: In-plane angular correlations of a particles (circles) and protons (triangles) 
measured in coincidence with the ER's Z = 16 (a) and 15 (b) in the 28 Si + 12 C reaction at 
E lab = 112 MeV in the angular range -115° < &{^ p < +115°). The proton correlations have 
been multiplied by a factor 10~ 3 for the sake of clarity. The arrow indicates the position 
of the IC detector at 0/ a b= -15°. On the abscissa, the positive angle refer to the opposite 
side of the beam from the direction of the ER detected in IC. The solid lines correspond to 
statistical-model calculations discussed in the text. 

Figure 7: In-plane angular correlations of a particles (circles) and protons (triangles) 
measured in coincidence with the ER's with Z = 17 (a), 16 (b), 15 (c), and 14 (d) in the 
28 Si + 12 C reaction at E lab = 180 MeV. in the angular range -115° < 0f a £ p < +115°). The 
proton correlations have been multiplied by a factor 10 -2 for the sake of clarity. The arrow 
indicates the position of the IC detector at Qi a b= -10°. On the abscissa, the positive angle 
refer to the opposite side of the beam from the direction of the ER detected in IC. The solid 
lines correspond to statistical-model calculations discussed in the text. 

Figure 8: Out-of-plane angular correlations of coincident a particles (circles) and protons 
(triangles) measured in the 28 Si + 12 C reaction at Ei ab = 180 MeV. The proton correlations 
have been multiplied by a factor 10~ 2 for the sake of clarity. The ER's with Z = 17 (a), 16 
(b), 15 (c), and 14 (d) are detected at Qi ab = -10°. The solid lines correspond to statistical- 
model calculations discussed in the text. 
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Figure 9: Energy-correlation plots between coincident a particles and S ER's produced 
in the 28 Si + 12 C reaction at E lab = 180 MeV. The heavy fragment is detected at <d s = -10° 
and the a-particle angle settings are given in the figure. The dashed lines correspond to 
different contours with their associated labellings discussed in the text. 

Figure 10: Experimental (a,c) and calculated (b,d) energy-correlation plots between 
coincident a particles and S ER's produced in the 28 Si + 12 C reaction at E iab = 180 MeV. 
The S is identified at Q s = -10° and the a particles are detected at 6 Q = +40°, and +70°, 
respectively. CACARIZO calculations are discussed in the text. 

Figure 11: Excitation-energy spectra calculated for the 28 Si + 12 C reaction at E/ ab = 180 
MeV for 8 Be (a) and 32 S in coincidence with the g.s.(b) and first excited level (c) of 8 Be. 
The solid line corresponds to the energy of the first excited state of 8 Be (3.08 MeV). The 
dashed lines correspond to an excitation energy in 32 S expected for an a-transfer process. 
The lower spectrum in (b) corresponds to the sum of the 4 individuals spectra. 
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